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Abstract 

5^ | In this paper we consider a Zi symmetrical 3-brane embedded in a 5-dimensional 

spacetime. We study the effective Einstein equation and acceleration condition in 
presence of the quintom dark energy fluid as the bulk matter field. It is shown 
that the time-dependent bulk quintom field induces a time-dependent cosmological 
qq ■ constant on the brane. In the framework of the DGP model, the effective Einstein 

equation is obtained in two different cases: i) where the quintom field is considered 
as the bulk matter field and the brane is empty and, ii) where the quintom dark 
energy is confined on the brane and the bulk is empty. We show that in both cases 

■ one could obtain a self-inflationary solution at late time in positive branch e = 1, 
OO ! and an asymptotically static universe in negative branch e = — 1. 
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1 Introduction 

Recent observations of type la supernova (SNIa) and WMAP [H [2] indicate that our uni- 
verse is currently undergoing an accelerating expansion, which confront the fundamental 
theories with great challenges and also make the researches on this problem a major en- 
deavor in modern astrophysics and cosmology. Missing energy density - with negative 
pressure - responsible for this expansion has been dubbed dark energy. Wide range of 
scenarios have been proposed to explain this acceleration while most of them can not 
explain all the features of universe or they have so many parameters that makes them 
difficult to fit. The models which have been discussed widely in literature are those which 
consider vacuum energy (cosmological constant) [3] as dark energy, introduce fifth ele- 
ments and dub it quintessence [I] or scenarios named phantom [5] with w < —1 , where 
w is parameter of state. A challenging issue is that the time-dependent dark energy gives 
a better fitting than a cosmological constant, and in particular the analysis of the prop- 
erties of dark energy from recent observations mildly favor models with w crossing -1 at 
redshift z « 0.2. Neither the quintessence nor the phantom alone can fulfill the transition 



* email: rezakord@ipm.ir 



1 



from w > — 1 to w < — 1 and vice versa. Although for k-essence [6J one can have both 
w > — 1 and w < — 1, it has been lately considered by Ref [7] that it is very difficult for 
k-essence to get w across —1 during evolving. But one can show [HI IS] that considering the 
combination of quintessence and phantom in a joint model, the transition can be fulfilled. 
This model, dubbed quintom, can produce a better fit to the data than more familiar 
models with w > — 1. In other words, the quintom model of dark energy represents a 
transition of dark energy equation of state from w > — 1 to w < — 1, or vice versa, namely 
from w < — 1 to w > — lis also one realization of quintom, as can be seen clearly in [ID] . 
Although the models with negative kinetic term are often plagued by instability, there 
are possibilities that these models might be phenomenologically viable if considered as 
effective field theories [TH [12] . 

An alternative way of explaining the observed acceleration of the late universe is to 
modify gravity at large scales. A well-studied model of modified gravity is the Dvali- 
Gabadadze-Porrati (DGP) braneworld model p3] where the brane is embedded in the 
flat bulk with infinite extra dimension. In this model gravity leaks of the 4-dimensional 
brane universe into 5-dimensional bulk spacetime at large scales. The inclusion of a 
graviton kinetic term on the brane recovers the usual gravitational force law scaling, 1/r 2 , 
at short distances, but at large distances it asymptotes to the 5-dimension scaling, 1/r 3 . 
Motivated by string/M theory, the AdS/CFT correspondence, and the hierarchy problem 
of particle physics, braneworld models were studied actively in recent years [H]-[I7J- In 
these models, our universe is realized as a boundary of a higher dimensional spacetime. 
The matter particles can not freely propagate in those large extra dimensions, but must be 
constrained to live on a 4-dimensional submanifold. The DGP model has a large scale/low 
energy effect of causing the expansion rate of the universe to accelerate. In almost all of 
works on braneworld models, the 5-dimensional bulk spacetime is assumed to be vacuum 
except for the presence of the cosmological constant, and the matter fields on the brane are 
regarded as responsible for the dynamics of the brane. However, from the unified theoretic 
point of view, the gravitational action is not necessarily the Einstein-Hilbert action. In 
fact, string theory tells us that the dimensionally reduced effective action includes not 
only higher-order curvature terms but also dilatonic gravitational scalar fields. Thus 
at the level of the low-energy 5-dimensional theory, it is naturally expected that there 
appears a dilaton-like scalar field in addition to the Einstein-Hilbert action [18] . Hence it 
is of interest to investigate how such a scalar field in the 5-dimensional theory affects the 
braneworld [13 [2D]. 

In this paper our main motivation is investigating the effects of the bulk quintom field 
on the evolution of the universe in the braneworld scenario and in the DGP model. We 
first review the braneworld scenario in presence of the bulk matter field in section 2. 
We study the acceleration condition for the universe with quintom dark energy in the 
bulk and show that the time-dependent bulk quintom field alters the brane as a time- 
dependent cosmological constant which is related to the quintom potential on the brane. 
In section 3 we obtain the generalized Einstein equation in the DGP model in presence 
of the tension and the bulk matter field. In the two following sections we investigate 
whether it is possible to have a late time accelerating phase on the brane when there is a 
quintom dark energy fluid in the bulk and the brane is empty; or inversely, when there is 
a quintom dark energy fluid on the brane and the bulk is empty. 
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2 Effective Einstein equation on the braneworld 



In this section we briefly review the braneworld scenario in presence of the tension and 
bulk matter field. In the braneworld scenario, our 4- dimensional world is described by 
a domain wall (brane) in 5-dimensional spacetime. We consider an ansatz for the 5- 
dimensional metric of the form 

ds 2 = -n 2 (t, y)dt 2 + a 2 (t, y)% j dx i dx j + b 2 (t, y)dy 2 , (1) 

where y is the coordinate of the fifth dimension and 7^ is a maximally symmetric 3- 
dimensional metric. We will use k to parameterize the spatial curvature and assume 
that the brane is a hypersurface defined by y — 0. We shall be interested in the model 
described by the action 

S = J d 5 xV^(^ 5) R -A + CT) + I <*V=g(-£ + C?*), (2) 

where OR is the scalar curvature of the 5-dimensional metric g AB , A is the bulk cosmo- 
logical constant, £ is the brane tension, k\ = SttG 5 , and q AB = 9ab ~ n A n B {n A is the unit 
vector normal to the brane and A,B = 0, 1, 2, 3, 5) is the induced metric on the 3-brane. 
The 5-dimensional Einstein equation can be written as 

{5) R AB - l -g AB ^R = k 2 (-Ag AB + T AB + S^ A 5 B 5(y b )), (3) 

here S(y b ) = T AB is the energy momentum tensor of the bulk matter and the last 
term corresponds to the matter content on the brane 

Sp, = -t,9nv + Tfa,. (4) 
The non-zero components of the 5-dimensional Einstein equation are 

3{-4-(- + b T ) + U- + -(- - t)) - - 2 } = ^ 2 (-a + t ° + s°5(y b )), (5) 

n z a a b z a a a a z 
1 „■ r a'a' n' . b'n' a\ a" n" , 

to 5 A- - + -) - rr(- + 2-) + - + — } + 

b z J a a n ba z n a a n 

+ 2-) - 2- + ^(-2- + ksi = k 2 5 (-A + t; + SjSfa)), (6) 

n z J a a n a b a n b J j j 

oX^ , a ' b ' a> 2 

3{ + -7 1 = h T 05, (7) 

n a a b a 

3{4(- + -) - --)+-)- 4> = ^- a + t ^ ^ 

ab z a n n z a a n a a z 

where primes indicate derivatives with respect to y, while dots derivatives with respect 
to t. Assuming a perfect fluid on the brane 

Ty = diag(-p b ,p b ,p b ,p b ), (9) 

and a quintom field in the bulk space containing the normal scalar field <j)(t, y) and negative 
kinetic scalar field a(t, y), with the Lagrangian expressed as the following form 

pnat = ^ g AB ^ _ ^ ^ + ^ (1Q) 



According to this action, the energy momentum tensor of the bulk quintom field is given 
by 

Tab = 4>,a4>,b - o,A(y,B - gAB(^g CD (4>,c(p,D - <t,cct,d) + V((f), a)). (11) 

In order to focus on the cosmological evolution on the brane we use the Gaussian normal 
coordinates (b(y,t) = 1) [21]. Thus the equations of motion of the scalar field <fi and a in 
bulk space are 



, )(p + n (— + 3— )<p +(p \-n — j- — = —t— < 

anna o<p d(p 



(12) 



(3 )<T + n ( h3— )a +a \+n = — 6{y b 

anna do da 



We are interested in studying the Einstein equation in presence of a quintom field in the 
bulk at the location of the brane. Without losing generality we choose n(t, 0) = 1 which 
can be achieved by scaling the time coordinate. As is well known, the presence of the 
brane leads to a singular term proportional to <5-function in y on the right-hand sides of 
the Einstein equations (jSJ) and ([6]) and the equation of motions (1121) . which have to be 
matched by singularity in the second derivatives in y on the left-hand side. Since all fields 
under consideration are symmetric under the orbifold symmetry Z 2 , these jumps in the 
first derivatives in y fix these first derivatives completely at y = 0. Here, these junction 
conditions read 

-U = -f(p, + 0, 

nVo = y(3^ + 2p 6 -0, (13) 



and 



\y=o 



1 5C^ at 



1 cfmat 

Using the components 00 and 55 of the Einstein equation in bulk space one can obtain 

F' = ^i(A-T >V"^T 5 Vd, (15) 

F = ^i(A-T 5 Va-^In 2 T 5 Va', (16) 
where F is a function of t and y defined by 



F(t, y) = - {a'af + ka 1 . (17) 



laa) , , 2 



Since the quintom field does not appear in the matter field Lagrangian on the brane 
(£™ ai ), Eq. (|14p implies that the quintom field is independent of y on the brane, namely 
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4>'\ y =o = &'\y=o — [19], [20]. Therefor the non- vanishing components of the quintom 
energy momentum tensor at the location of the brane are 

Tl=p B = l -4> 2 - l -a 2 -V{^a), (18) 

Tq5 = #' - = 0, 

As one can see, in the case of y independent bulk quintom field T05 vanishes. It means 
that there is no flow of matter along the fifth dimension. Using ( [18]) one can solve Eq.( fT5l) 
which leads to the first integral of the 00 component of Einstein equation as 

o a 4 a z n z cr ar 

where C is a constant of integration which is usually referred to dark radiation [22] . 
Substituting the junction conditions (jT3l) into above equation, we arrive at the generalized 
Friedmann equation on the brane as 

h h 2 h 2 k 4 h 2 h 4 C 

h 2 + 4 = i ( A + i a + + ip* + 1^ 2 + ^ ( 2 °) 

a u o 18 o do a 4 



here = - is the Hubble parameter. As one can see from Eq. ff20|) . in the absence of the 
bulk matter field, the cosmological constant and the brane tension, the equation gives rise 
to a Friedmann equation of the form H oc pb instead of H oc ^/pb which is inconsistent 
with cosmological observation. This problem can be solved by either considering the 
cosmological constat and tension on the brane or considering a matter field in the bulk 
[23| [2l]. Recalling the junction conditions (|T3l) . the 05 component of Einstein equation 
and field equations ([12]) on the brane take the following form respectively 

p b + 3H(p b + p b ) = 0, (21) 
+ 3-0+^^ = 0, (22) 

z +3 * & _ mM = , (23) 

a acr 

It should be noted that if scalar field and a satisfy the field equations ([22]) and ( [23]) 
respectively, the bulk energy momentum tensor is automatically conserved and we have 

p B + 3H{p B +p B ) = 0. (24) 

We are interested in studying the acceleration condition for a universe with the quintom 
field in the bulk. The condition for acceleration can be obtained from ( 120]) by using the 
conservation equation of the brane and bulk matter field ([2T]) and ([24]) 



\ = f (A + §£ 2 ) - + 3p fc ) - §( PB + 3 PB ) - M (2p g + 3 P ^) - £ (25) 
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To study the role of the bulk quintom field in the late time acceleration phase on the brane, 
we ignore the effect of tension, brane matter, cosmological constant and dark radiation 

\ = -%{Pb + *Pb). (26) 

Thus the acceleration condition for a universe with quintom dark energy in bulk is 

PB<-y, or 4> 2 -& 2 <V(cf>,a). (27) 

Now we consider the 55 component of Einstein equation at the position of the brane which 
leads to the Raychaudhuri equation 

^ + ^ 2 + ^ = |(A+ f£ 2 ) - |(£(3p 6 - Pb ) + Pb (3p b + Pb )) - Mi* (28) 

Using Eq. ([20]) one can rewrite the above equation as 

h P P P P P C P 

5 = |(A + h e) _ k, {pb + 3pi ) - h„ B - h (2pi + 3ptPb) - | - (29) 

Comparing Eq. (l25l) with (|29|) provide a constraint on the bulk energy momentum tensor 

(3p B - p B ) = 4T 5 5 , (30) 

which for the quintom field with the energy momentum tensor ([18]) leads to 

2 -a 2 = O. (31) 

It means that the time-dependent bulk quintom field influences the brane like a time- 
dependent cosmological constant which can be written in terms of the quintom potential 
energy, pb = —ps = — V{4>) [2S]- For a particular solution of (j3Tl) in which <p and a are 
both constant on the brane, we arrive at the natural cosmological constant induced by 
the time-dependent bulk quintom field. 



3 Generalized Einstein equation in the DGP braneworld 

In the DGP model, which provide a simple mechanism to modify gravity at large dis- 
tances, it is supposed that a 3-dimensional brane is embedded in a flat 5-dimensional 
bulk. This model predicts that 4-dimensional Einstein gravity is a short-distance phe- 
nomenon with deviations showing up at large distances. The transition between four- and 
higher-dimensional gravitational potentials in the DGP model arises as a consequence of 
the presence of both brane and bulk Einstein terms in the action. The DGP model in- 
cludes a length scale below which the potential has usual Newtonian form and above 
which the gravity becomes 5-dimensional. The cross over scale between the 4-dimensional 

k 2 o 

and 5-dimensional gravity is r c = ^\ in which /i = 87^4. In this framework, existence 
of a higher dimensional embedding space allows for the existence of bulk or brane matter 
which can certainly influence the cosmological evolution on the brane. Now we proceed to 
obtain the generalized DGP model in which both bulk cosmological constant A and brane 
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tension £ are non-zero. We consider the model described by the gravitational bulk-brane 
action 

S = J d 5 x^(^R -A + CT l ) + J d'x^i^R -i + £™*), (32) 

where ^'R is the Ricci scalar of the induced metric q^ u . The 5-dimensional Einstein 
equation takes the form 

{5) Rab - = k%(-Ag A B + T AB + S^ A 5 B 5(y b )). (33) 

Here Tab is the energy momentum tensor of the bulk matter, A is the cosmological 
constant of the bulk spacetime and the energy momentum tensor on the brane is given 
by 

S^u = ~Cq^ + V ~ ^~ 2 U^u, (34) 

the last term is the contribution coming from the scalar curvature of the brane with the 
non-vanishing components given by 



U 00 = 3(# 2 + W 



J7« = (^(-H 2 + 2H^-2^)-k) lij . 
n z n a 



(35) 



From 00 and ij components of the Einstein equation (133]) we find the following junction 
conditions which simply relate the jumps of derivatives of the metric across the brane to 
the stress tensor inside the brane 

-\ y=0 = -^(p 6 + + r c (H 2 + 4), 
a o a z 

n'\y=o = ^-(Sp b + 2p b - + r c (-H 2 + 2--\). (36) 
6 a a 1 



In the following two sections we study the Einstein equation ( 13311 when the quintom field 
is considered as the bulk matter field, and a quintom dark energy confined on the brane, 
respectively. 



3.1 Quintom field in bulk space 

We consider the quintom field as the bulk matter field with Lagrangian expression in 
Eq. (IT0l . Integrating the equation of 00 component of fl33|) around y = and using 
junction conditions fl36l) . we arrive at the generalized (first) Friedmann equation 

(l + |^)(^ 2 +4)-|(A + |a-|^-|p B +-^ = ^(Vp,+3(^ 2 +4)) 2 . (37) 
o/i ar 6 6 18 6 a 4 36/i 4 a z 

The brane-Friedmann equation ( l20l) can be derived from above equation by letting /x go to 
infinity. We are interested to study the effect of the quintom field on the brane. Ignoring 
the cosmological constant A, the brane tension and the matter field on the brane, Eq. (1371) 
can straightforwardly be rewritten as 
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here ps is the energy density of the bulk quintom field on the brane derived in Eq.flTSJ). 
The two different possible e namely e = ±1, correspond to two different embeddings of 
the brane into the bulk spacetime [26]. Since the bulk quintom field satisfies the usual 
energy momentum conservation law on the brane (|24|) . we have ps oc a~ 3 ( w+1 ) [to is the 
state parameter). Integrating Eq. (1381) for k = and uo > — 1 where <fi > a, shows that the 
scale factor a diverges at late timdjj (See Figure 1). Thus the energy density of the bulk 
matter goes to zero at late time and reaches a regime where it is small in comparison with 
1/r 2 . In the case u < — 1 where <f) < &, integrating Eq. (l38l indicates a vanishing scale 
factor a at late time, so the matter density goes to zero and we could use the assumption 
k 2 ps <C Therefore, in the DGP model with a quintom dark energy fluid in the bulk 

space, one can expand the Einstein equation ([3"51) under the condition that k\pB <C 1/r 2 
for all range of u. At zero order and for spatially flat metric, two different results depend 
on the value of e can be derived. Considering the case e = — 1 yields 

H 2 = 0, (39) 

which describes an asymptotically static universe. In the other case we take e = 1 which 
leads to 

H 2 = ^r, or ait) ocexpf— ). (40) 

This provides the self-inflationary solution at late time which is the most important aspect 
of the DGP model. Therefor, the late time behavior of the universe does not alter even if 
we ignore the matter field on the brane and consider a model of the universe filled with 
the bulk quintom dark energy. 

(a) (b) (c) (d) 

180,000 -j 
160,000 - 
140,000 - 
120,000 - 
100,000 - 
80,000 - 
60,000 - 
40,000 - 
20,000 - 

<H \ 

-100 -50 50 100 -50 50 100 -100 -50 50 100 -50 50 100 

a a a a 

Figure 1: These figures show the evolution of / as a function of scale factor. 



1 The integration of Eq. (|38jl leads to I = J — , a = -J^- J dt in which m > when to > — 1 

2fc 2 r 2 

and to < when u> < — 1 and p = — The variation of I against a in different value of m and e are 
plotted in Fig.l 
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3.2 Quintom field on the brane 



In this section we ignore the bulk matter and consider a quintom dark energy confined 
on the brane in the DGP model with the Lagrangian expression as 

C rnat = 1^(0 ^ _ a ^ y ) + V(<j>, a). (41) 

The energy momentum tensor of the quintom field on the brane is given by 

t^u = 4>,n4>,v - o>°> - ^(^ Q/3 (</>,a</\/3 - v,^,?) + v(4>, a)). (42) 

In absence of the bulk matter field and the brane tension, the generalized Friedmann 
equation (1371) leads to [27] 

^4=^W i+ ¥^ 2) - (43) 

in which is quintom energy density obtained by ( ]42l . Since the energy momentum 
tensor on the brane is conserved, we could apply the strategy used in the previous section 
to study the late time cosmology on the brane. Integrating Eq. (|43|) for a spatially flat 
spacetime yields the same result as the previous section. It indicates that for uj > — 1 
the scale factor a diverges at late time [27J and for uj < —1 the scale factor a vanishes 
at late tinted (Figure 2). Thus the energy density of quintom dark energy goes to zero 
for late time and reaches a regime where it is small in comparison with \jr 2 c . Expanding 
the equation fT43|) under the condition /i 2 p& <C 1/r 2 provides an asymptotically static 
universe, H — 0, in the case e = — 1 and a self-accelerating phase, H — — , in the case 
e = 1. Therefore the presence of quintom dark energy on the brane or in the bulk dose 
not change the late time behave of the universe. In both cases for all range of uj {uj < — 1 
or uj > —1), one can derive the self- accelerating universe at late time. 



(a) (b) (c) (d) 




Figure 2: These figures show the evolution of i 7 as a function of scale factor. 



4 Conclusion 

We have studied the cosmology of a Z 2 symmetrical 3-brane embedded in a 5-dimensional 
spacetime including a quintom dark energy fluid in bulk space. In the braneworld sce- 
nario, we derived the acceleration condition due to the bulk quintom field. It was indi- 
cated that the time-dependent bulk quintom field alters the dynamics on the brane as 

2 The integration of Eq.(H51) leads to /' = f d ? = ^- f dt in which m > when uj > -1 and 

J o(l+«Vl+Jr) lr " J 

m < when uj < — 1 and p = 3 c , 
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a time-dependent cosmological constant which can be derived in terms of the quintom 
field potential. It means that to have an accelerated expansion phase on the brane, the 
potential energy due to the bulk quintom field must be a positive function of time. 

In the DGP model, when an intrinsic curvature term is added on the brane, we have 
obtained the generalized Einstein equation, where both bulk and brane matter field are 
non-zero. Cosmology on 3-brane have been studied in two different cases. In first case, we 
have considered the quintom field as the bulk matter field and ignored the brane matter. 
It was shown that two different solutions are obtained for two different embeddings of the 
brane. In negative branch where e = — 1, the generalized Friedmann equation describes an 
asymptotically static universe, and in positive branch where e = 1, we obtained H — — , 
which provides a self-inflationary solution at late time. In second case, we have considered 
a quintom dark energy confined on the brane in the DGP model, and studied the solution 
of the generalized Friedmann equation, when the bulk matter field and the brane tension 
were ignored. Similar to the first case, here also two different results can occur, H = 
for e = — 1 and a self-inflationary solution at late time, H = — , for e = 1. Therefor, the 
late time behavior of the universe does not alter even if we ignore the matter field on the 
brane and consider a model of the universe filled with the bulk quintom dark energy, or 
vice versa, we ignore the bulk matter field and consider only the quintom dark energy on 
the brane. 

Finally we should stress on the ghost instabilities present in the self-accelerating branch 
of this DGP-inspired model. The self-accelerating branch of the DGP model contains a 
ghost at the linearized level [28]. Since the ghost carries negative energy density, it leads 
to the instability of the spacetime. The presence of the ghost can be attributed to the 
infinite volume of the extra-dimension in DGP setup. When there are ghosts instabilities 
in self-accelerating branch, it is natural to ask what are the results of solutions decay. As 
a possible answer we can state that since the normal branch solutions are ghost-free, one 
can think that the self-accelerating solutions may decay into the normal branch solutions. 
In fact for a given brane tension, the Hubble parameter in the self-accelerating universe is 
larger than that of the normal branch solutions. Then it is possible to have nucleation of 
bubbles of the normal branch in the environment of the self- accelerating branch solution. 
This is similar to the false vacuum decay in de Sitter space. However, there are arguments 
against this kind of reasoning which suggest that the self-accelerating branch does not 
decay into the normal branch by forming normal branch bubbles ( see [28] for more 
details). It was also shown that the introduction of Gauss-Bonnet term in the bulk does 
not help to overcome this problem [25]. In fact, it is still unclear what is the end state 
of the ghost instability in self-accelerated branch of DGP inspired setups. On the other 
hand, quintom scalar fields and induced gravity in our setup provides a new degree of 
freedom which requires special fine tuning and this may provide a suitable basis to treat 
ghost instability. It seems that in our model this additional degree of freedom has the 
capability to provide the background for a more reliable solution to ghost instability due 
to wider parameter space. 

5 Acknowledgment 

We would like to thank an anonymous referee for helpful comments on the paper. 



10 



References 



[1] S. Perlmutter et al, Nature (London), 391, 51, (1998); Knop. R et al., Astroph. J., 
598, 102 (2003); A. G. Riess et al., Astrophy. J., 607, 665(2004); H. Jassal, J. Bagla 
and T. Padmanabhan, Phys. Rev. D, 72, 103503 (2005). 

[2] D. N. Spergel et al. [WMAP Collaboration], ApJS, 170, 377 (2007); L. Page et al. 
[WMAP Collaboration], ApJS, 170, 335 (2007); G Hinshaw et al. [WMAP Collabo- 
ration], ApJS, 170, 288 (2007); N. Jarosik et al. [WMAP Collaboration], ApJS, 170, 
263 (2007). 

[3] For review on cosmological constant problem: P. J. E. Peebles, B. Ratra, Rev. Mod. 
Phys., 75, 559-606, (2003); J. Kratochvil, A. Linde, E. V. Linder, M. Shmakova, 
JCAP, 0407, 001, (2004). 

[4] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett., 80, (1998) 1582; I. 
Zlater, L. Wang and P. J. Steinhardt, Phys. Rev. Lett., 82, (1999), 896; M. S. Turner 
and M. White Phys. Rev. D, 56, (1997), 4439. 

[5] R. R. Caldwell, Phys. Lett. B 545, 23, (2002); R. R. Caldwell, M. Kamionkowsky and 
N. N. Weinberg, Phys, Rev, Lett, 91, 071301, (2003); S. Nojiri and S. D. Odintsov, 
Phys. Lett., B 562, (2003), 147; S. Nojiri and S. D. Odintsov, Phys. Lett., B 565, 
(2003), 1; M. R. Setare, Eur. Phys. J. C 50, 991, (2007). 

[6] C. Armendariz-Picon, V. Mukhanov and P. J. Steinhardt, Phys. Rev. Lett. 85, 4438 
(2000); Phys. Rev. D 63, 103510 (2001); T. Chiba, T. Okabe and M. Yamaguchi, 
Phys. Rev. D 62 (2000) 023511. 

[7] A. Vikman, Phys. Rev. D 71, 023515 (2005). 

[8] B. Feng, [isTro-ph/060215"6 



[9] Z. K. Guo et al, astro-ph/0410654 



[10] Gong-Bo Zhao, Jun-Qing Xia , Bo Feng , Xinmin Zhang , astro- ph/0603621; Jun- 
Qing Xia, Gong-Bo Zhao, Bo Feng , Xinmin Zhang , astro-ph/0603393; Bo Feng, 
Mingzhe Li , Yun-Song Piao , Xinmin Zhang, Phys. Lett. B634, 101, (2006); M. R. 
Setare, Phys. Lett. B 641, 130 (2006); J. Sadeghi, M. R. Setare, A. Banijamali, F. 
Milani, Phys. Lett. B662, 92, (2008); M. R. Setare, J. Sadeghi, A. R. Amani, Phys. 
Lett. B660, 299, (2008); M. R. Setare, E. N. Saridakis. la7Xiv:0802.2595l [hep-th], 

[11] S. M. Carroll, M. Hoffman and M. Trodden, Phys. Rev. D 68, 023509 (2003). 

[12] Gong-Bo Zhao, Jun-Qing Xia , Mingzhe Li , Bo Feng, Xinmin Zhang , Phys. Rev. 
D72, 123515, (2005). 

[13] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000). 

[14] P. Horava and E. Witten, Nucl. Phys. B460, 460 (1996). 

[15] L. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 



11 



[16] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B436, 
257 (1998); N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B436, 263 
(1998). 

[17] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); L. Randall and R. 
Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 

[18] A. Lukas, B.A. Ovrut, K.S. Stelle and D. Waldram, Phys. Rev. D59, 086001 (1999); 
Nucl. Phys. B552, 246 (1999); Y. Himemoto, and M. Sasaki, Phys. Rev. D63, 044015 
(2001). 

[19] U. Ellwanger, Eur.Phys. J.C25, 157 (2002). 

[20] Sean M. Carroll and L. Mersini, Phys. Rev. D64, 124008, (2001). 

[21] P. S. Apostolopoulos and N. Tetradis, Class. Quant. Grav. 21, 4781 (2004) 
[arXiv : hep-th/ 0404 1 5 1 ; P. S. Apostolopoulos, N. Brouzakis, E. N. Saridakis and 
N. Tetradis, Phys. Rev. D 72, 044013 (2005) [arXiv:hep-th/0502115] . 

[22] R. Maartens, Phys. Rev. D 62 (2000) 084023 

[23] J. M. Cline, C. Grojean and G. Servant, Phy. Rev. Lett. 83, 4245 (1999). 
[24] P. Kaonti, I. Kogan, K. Olive and M. Pospelov, Phys.Lett. B 468, 31(1999). 
[25] C. Wetterich, Astron.Astrophys.301, 321,(1995). 

[26] M. Cvetic, S. Grimes and H. H. Soleng, Phys. Rev. D48 (1993) 2613 ; G. W. Gibbons, 
Nucl. Phys. B394 (1993) 3. 



[27] C. Deffayet, Phys. Lett. B 502 (2001) 199 |hep-th/0010l86] . 



[28] K. Koyama, Class. Quantum Grav. 24, R231 (2007) [arXiv:hep-th/0709.2399]. 
[29] C. de Rham and A. J. Tolley, J CAP 0607, 004 (2006) |arXiv:hep-th/0605122| . 



12 



